Fragile-glass behavior of a short range p-spin model 
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In this paper we propose a short range generalization of the p-spin interaction spin-glass model. The model is well 
suited to test the idea that an entropy collapse is at the bottom-line of the dynamical singularity encountered in 
structural glasses. The model is studied in three dimensions through Monte Carlo simulations, which put in evidence 
fragile glass behavior with stretched exponential relaxation and super- Arrhenius behavior of the relaxation time. Our 
data are in favor of a Vogel-Fulcher behavior of the relaxation time, related to an entropy collapse at the Kauzmann 
temperature. We however encounter difficulties analogous to those found in experimental systems when extrapolating 
thermodynamical data at low temperatures. We study the spin glass susceptibility investigating the behavior of the 
correlation length in the system. We find that the the increase of the relaxation time is not accompanied by any 
growth of the correlation length. We discuss the scaling properties of off-equilibrium dynamics in the glassy regime, 
finding qualitative agreement with the mean-field theory. 

72.10.Nr, 64.60. Cn 



I. INTRODUCTION 

The glassy state is very common in natureS. When it is 
reached from the liquid phase lowering the temperature, 
one finds a dramatic increase of the relaxation time, and 
off-equilibrium phenomena can not be avoided on exper- 
imental time scales. This leads to non analytic behavior 
of the thermodynamic quantities, with a "transition tem- 
perature" that depends on the cooling rate. Despite its 
ubiquity, the basic mechanisms underlying the common 
features, as well as the peculiarities of the glassy behavior 
in different systems are yet to be clarified. 

One of the most suggestive ideas in the glass theory, 
proposed long time ago by Gibbs and Di Marziocl, relates 
the increase of the relaxation time, and the observed fi- 
nite time singularities to the existence of a thermody- 
namic transition at the Kauzmann temperature where 
the configurational entropy collapses to zercQ. Soon af- 
ter, in a refinement of the argument, Adams and Gibbs 
argued in favor of a Vogel-Fulcher singularity in the re- 
laxation time. 

Disordered systems have been proposed as paradig- 
matic models in which glassy phenomena can be stud- 
ied in a nutshell, and theoretical ideas tested on micro- 
scopic modelsy. This is due to the fact that in disordered 
systems, the glassy state already appears in mean-field 
theory. The natural separation of the variables among 
"quenched" and "annealed" allow for the successful use of 
powerful techniques as the replica method for staticB and 



functional methods in dynamicsi. In fact, a satisfactory 
mean-field theory of disordered systems for the statical as 
well as Jpr the equilibriums and off-equilibrium dynamics 
existaO'Q. Recently, examples of mean-field deterministic 
models with glassy behavior very similax to the one of 
disordered systems have been displayed^. This points 
in the direction that common mechanisms could lead to 
the glassy behavior of disordered and non-disordered sys- 
tems. 

The simplest example in which the Gibbs-Di Marzio 
collapse— .occurs is the Random Energy Model of 
Derridali9, and is a common feature to all systems with 
a "discontinuous glassy transition" , or technically "one- 
step replica symmetry breaking" , where the Edwards- 
Anderson parameter!!!! undergoes a discontinuity!]. Ex- 
amples of such models., are the Potts Glasal3, the p- 
spin interaction modeled and a model of manifolds in— a 
disordered media with short-range correlated disorder!!.! 
This class of systems has been proposed by Kirkpatrick, 
Thirumalai and Wolynes as simple toy models for the 
structural glass transition!!!. Notably, the study of the 
Langevin dynamics of the spherical version of these mod- 
els shows that there the Mode Coupling TheorytLa is ex- 
act, and displays a dynamical singularity of kind B in 
the Gotze classification!!^!. In fact, recent progresses in 
the comprehension of the dynamics of mean-field disor- 
dered systemsQIa allowed for an extension of the-Mode 
Coupling Theory to the-broken ergodicity phaselljtl. 



Many studies 



have pointed out the exis- 
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tence of a temperature To where, despite the fact that 
no singularity is observed in the free energy, there is 
a statical breaking of ergodicity into an exponentially 
large number of metastable states. A thermodynamic 
singularity is present at a temperature Tq smaller then 
To- As a genuine mean- field theory, the Mode Coupling 
Theory neglects the activated (droplet) processes that 
in finite-dimensional systems are responsible for the de- 
cay of metastable states in a finite time. Kirkpatrick 
and Wolynes have recently stressed how the inclusion of 
these processes can restore ergodicity for T c < T < T D , 
and give rise to a generalized Vogel-Fulcher singularity 
at the static transition temperature T(ja. The argument 
has been confirmed and refined by Parisi with a theo- 
retical .calculation based on the potential theory in spin 
glassesE3. 

The aim of this paper is to test this idea in a finite- 
dimensional disordered model where metastable states 
should be present, but destabilized by activated pro- 
cesses. In section 2 we propose a finite dimensional analo- 
gous of the p-spin interaction model in the case p = 4. We 
believe that, as usual, the mean-field limit is recovered for 
high dimensionality. In section 3 we study the thermo- 
dynamics in the high-temperature regime through Monte 
Carlo simulations, that demonstrate that the model be- 
haves as a fragile glass. 

A second aspect of our work concerns the off- 
equilibrium dynamics deep in the glassy phase. There, 
the properties of the system depend on the thermal his- 
tory, and time translation invariance does not holdEl 
The off-equilibrium Mode Coupling Theory predicts scal- 
ing relations and a definite pattern of violation of the 
Kubo fluctuation-dissipation relation. In section 4 we 
study the dynamics in this regime, showing the consis- 
tency of the aforementioned scenario. Finally, the con- 
clusions are drawn. 



II. THE MODEL 

The p— spin modelEHil is defined by the long-range 
Hamiltonian 
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where the couplings Jj 1; 



are independent Gaus- 



sian variables with zero mean, and variance J? „• „ = 
p\/(2N p - 1 ). The spins a u i = 1,...,N can be taken 
as Ising variables, or as real variables subjected to the 
spherical constraint °f = A0 The case p = 2, 

that corresponds to the Sherrington Kirkpatrick modeL 
has a glassy transition with continuous order parameters 
and will not be considered in this paper. For p > 3 both 
in the Ising and in the spherical case the transition is 
discontinuous, and the properties of the model are the 
ones of interest in this paper. As a finite dimensional 



analogous of the model (|]) in the case p = 4 we take a 
spin system with interacting Ising spins <7j arranged on 
the sites of a D dimensional square lattice with periodic 
boundary conditions. The Hamiltonian is defined as 
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where the sums runs over all the plaquettes □ of 
the lattice. Each spin belongs to 4 ( ^ ^ plaquettes. 

Each plaquette □ gives a contribution —Jn Ilien °» ~ 
Jn (7° o"2 CT3 CT4 , where the variables Jn are chosen as in- 
dependent Gaussian variables with zero mean and unit 
variance. Note that, in generic dimension, the model is 
not invariant under Z 2 gauge transformation as it would 
be if the spins were located on the links (instead that 
on the vertices) of the plaquettes. Thanks to this, there 
is no Elitzui theorem preventing non-zero global order 
parameterstj. In fact there .are no local symmetry op- 
erations leaving H invariants. However the Hamilto- 
nian is invariant under the contemporary inversion of all 
the spins that belong to any hyperplane of dimension 
D — 1 orthogonal to one of the Cartesian axes. It is easy 
to check that these operations do not change the sign 
of any of the plaquettes. The degeneracy due to this 
symmetry (plane inversion symmetry in the following) is 
2-DL-(D-i)^ an( j can kg removec [ e g fixing the spins on 

the Cartesian axes to arbitrary values. This exponen- 
tially large degeneracy of the states is also present in the 
ferromagnetic version of the model (Jn = 1). We think 
that this could lead to a very interesting spinodal dynam- 
ics, above the lower critical dimension D = 2. We con- 
centrate here to the disordered model, leaving the study 
of the ferromagnetic case for future work. A ferromag- 
netic model with four spin-interactions at the vertices of 
plaquettes was studied inc3 in connection with random 
surfaces physics. In that case also pair interactions that 
removed the plane inversion symmatry where present in 
the Hamiltonian. 

In the limit of infinite dimension, where the number 
of plaquettes to which a spin belongs tends to infinity, 
one can expect that, modulo the symmetry, the model 
is equivalent to the p— spin model for p = 4. Models 
with different p's could be easily constructed for other 
lattices, e.g.pthe case p — 3 would correspond to a trian- 
gular latticecll. 

It is worth at this point to present a brief qualitative 
review of the result^pf Ihc mean-field theory, based on 
the Hamiltonian (0)q1!J'EI The study of the thermody- 
namics of this system leads to the following results. At 
high temperatures the system is paramagnetic and er- 
godic. At a temperature To the ergodicity breaks down, 
and an exponentially large number (exp(7V£(T))) of pure 
states (ergodic components) separated byJparriers of or- 
der N contribute to the partition functional. This transi- 
tion occurs without singularities in the free energy, which 
is equal to the free-energy per state F in , plus an entropic 
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contribution — TH(T) coming from the multiplicity. The 
quantity S, the configurational entropy, is a decreas- 
ing function of the temperature, and at a temperature 
Tq < Td vanishes. At Tq there is a thermodynamic 
phase transition with singularities in the free-energy. 

On the other hand, the study of the off-equilibrium 
dynamics after a sudden quench from a high tempera- 
ture shows that the large-time limit of various dynamical 
quantities is non analytic and the singularity in the dy- 
namics at Td is suppressed. The characteristic time scale 
r for these processes to restore ergodicity has been re- 
cently estimated in Potts-glass models using heuristic ar- 
guments by Kirkpatrick and WolynesH, and substantiated 
using a droplet argument in replica space by ParisiE3. 
They find a generalized Adam-Gibbs relation of the kind 



exp 



C 
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where C is a constant and 7 = D — 1. As the configu- 
rational entropy vanishes linearly near 

Tc, eq-(|) results 

in a generalized Vogel-Fulcher law r ~ exp ( ^ T J^ c y ^ ■ 

The value 7 = 2 in D = 3 is at variance with the usual 
value 7 = 1 used to fit the experimantal data. However, 
in the case of the present model the value of 7 = D — 1 
in (|^) should be lowered due to the plane inversion sym- 
metry. We do not know if this would result in 7 = D — 2, 
and we leave the investigation of this point for future 
work. As the matter of fact, simple (and trivial) results 
are obtained for the statics in D = 2. In that case, 
one can show that in the high temperature expansion 
only diagrams involving a number of spins proportional 
to L or higher, and hence irrelevant for L — > 00, are 
present. Accordingly the free-energy per spin is found 
F = -Tf -^exp(-J 2 /2)log(2cosh(/3J)). We will see 
that the relaxation time follows a simple Arrhenius law 
in this case. The lowest dimension at which one can ex- 
pect non trivial thermodynamical results is D = 3, where 
one can see that frustration is present. The study of the 
properties of the three-dimensional model through Monte 
Carlo simulations and the comparison with the results of 
the theory, will be the subjects of the rest of this paper. 
Some results for the two dimensional case will also be 
mentioned. 



III. THERMODYNAMICS 

In order to investigate the questions posed in the previ- 
ous sections, we have performed Monte Carlo simulations 
of the model in three dimensions, using a standard serial 
single spin-flip heat-bath algorithm. The signature of 
glassy behaviour is easily seen in simulations of cooling 
experiments. In figure 1 we plot the energy as a function 
of the temperature for different cooling rates. We clearly 
see a change of behaviour corresponding to a jump of the 
specific heat around T ~ 0.7, where the systems fails to 



reach equilibrium within the observation time. The in- 
set shows that the "transition temperature" as well as the 
value of the energy at which the system freezes are depen- 
dent on the cooling rate. Since we expect the equilibrium 
entropy to be the relevant quantity for the transition, wc 
integrated the high-temperature energy data in /3 to get 
the entropy, taking into account data from (3 = 0.01 up 
to = 1. 

In figure 2 we present the results of this operation, 
togheter with some rational function fits of the data 
points. The functional form that we have chosen to fit 
that data are: 



log(2) + a/3 2 
1 + 6/3 2 
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log(2) + a/3 2 
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Even functions of [3 have been chosen coherently with the 
fact that the high temperature expansion of S only con- 
tains even powers of [3. The fits are roughly of the same 
quality at high temperature and are both consistent with 
the entropy collapse scenario. Extrapolating at low tem- 
perature, one can estimate the point where the entropy 
vanishes to be in the range 0.3 < T c < 0.6. This range 
of temperatures should not be considered as more than a 
mere indication of what it could happen. First of all, the 
extrapolation at low temperatures by means of the func- 
tions (^) is highly arbitrary. Second, we expect only the 
configurational entropy, and not the total entropy (which 
take an intra-state contribution), to become negative at 
the transition. The impossibility to disentangle the two 
contributions adds incertitude to the critical tempera- 
ture estimate. However, both the mean-field theory and 
the off-equilibrium dynamics of section 4 suggest that 
in all the low-temperature region the Edwards- Anderson 
parameter is close to 1, correspondingly the intra-state 
entropy is rather small, and possibly much smaller than 
the configurational one not too close to T c . 

Coherent informations are obtained from the study 
of the equilibrium dynamics at high temperatures. 
The analysis of the autocorrelation function C(i) =< 
a(t)a(0) > shows how the relaxation follows a stretched 
exponential law of the type, 



C(t) = exp(-(t/r) b ) 



(5) 



As shown in figure 3 our data are pretty well fitted by 
the previous relaxation law all along the relaxation. 

From that fit we are able to extract the tempera- 
ture dependent relaxation time r(T) and the exponent 
b(T). The behavior of these parameters as a function 
of the temperature is depicted in figure 4. The relax- 
ation time is consistent with a Vogel-Fulcher law of the 
kind t = Aexp(B/(T - T )) with T ~ 0.58, A = 0.5, 
B = 10 (continuous line in figure 4). But, analogously 
with what happens with experimental data, it is also 
well fitted by a zero temperature singularity of the type 
t = Aexp(B/T 2 ) with A ~ 11.5, B ~ 0.69. A four pa- 
rameter fit of the form r = Aexp(B/(T - T ) 7 ) fits the 
data with 7 = 1.01, while posing by hand 7 = 2 as in 
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(||) leads to Tq = 0.001. The indication of a thermody- 
namic transition by the divergence of the relaxation time 
at To is supported by the fact that it is in the range of 
temperatures where the extrapolated entropy vanishes. 
The exponent b is also depicted in the inset of figure 4. 
It appears to be linear with the temperature and of or- 
der .5 in the low temperature phase. Unfortunately we 
have no evidence stronger than this in favour of the finite 
temperature singularity. The best we can say is that our 
numerical data are consistent with the glass transition 
scenario as much as laboratory experiments on glasses 
give support to this singularity (with the difference that 
laboratory experiments can explore a larger window of 
time than in the numerical experiments). However the 
picture we get is coherent with the theoretical relation 
between the relaxation time and the configurational en- 
tropy. We have related the two quantities as suggested by 
eq. (|^) . In figure 5 we plot the logarithm of the relaxation 
time versus the inverse entropy, and we see that the data 
fall quite well on a straight line indicating the validity of 
(|J) with 7 = 1. This linear relation is again consitent 
with a small value of the intra-state entropy, that could 
be estimated of the order of Soo = 0.04, value obtained 
from the extrapolation of the fit in fig. 5 to r — > oo. 

It is worth at this point to mention the results of an 
analogous analysis for the two dimensional system, where 
the thermodynamics is trivial. Simulations performed 
in that case indicate that also in D = 2 there exists a 
low-temperature regime where the equilibrium correla- 
tion function behaves as a stretched exponential (we ob- 
serve that stretched exponential behavior is also observed 
in the one-dimensional Ising model at low temperatures 
in the intermediate time regimeEa). However, the relax- 
ation time grows much slower at low temperature than 
in the three dimensional case. In fact, very good fits are 
obtained by simple Arrhenius forms r ~ exp(B/T), i.e. 
the minimal increase to be expected in any system with 
discrete spin variables. 

We note that the two-step relaxation characteristic of 
many structural glasses is not seen here. 

To clarify further the picture about the thermodynam- 
ics of the three dimensional model, we have investigated 
the possibility of a growing static correlation length. This 
is done in disordered spia. system through the analysis of 
the correlation functioned 



G(x) = (a T a x T x ) - exp(--) (6) 

among two different replicas <Tj and r, with same Jo's, 
and evolving with independent thermal noise. To avoid 
problems with the symmetry of the system wc have fixed 
all the spins <Ji and Tj belonging to the Cartesian axes. 
In the thermodynamic limit this should be an effective 
way to eliminate the degeneracy. 

In our simulations we have not found any evidence of 
a growing correlation length, which appears to be inde- 
pendent of temperature and smaller than 1. We conclude 



that the divergence of the relaxation time is not accom- 
panied by a divergent correlation length. 

An alternative way to explore the existence of a di- 
vergent correlation length is to compute directly the in- 
tegral of the correlation function G(x) over the whole, 
space. This quantity yields the spin-glass susceptibilityEll 
defined by, 

XSG = N{W)-{V)f) (7) 

where N = L 3 is the size of the system. For vanishing 
correlation length this quantity is equal to 1 and tends to 
increase whenever theije,is spatial ordering and the corre- 
lation length increases^. We have simulated several sys- 
tem sizes L = 4, 5, 6, 7 from T = 3.0 down to T = 0. Data 
is shown in figure 6. As one can expect, there are serious 
thcrmalization problems especially at low temperatures. 
Nevertheless, it clearly emerges form these results that, 
in the region of temperatures where we are in equilib- 
rium (for instance, above T = 1.2 where the equilibrium 
relaxation time is smaller than the thermalization time 
we made the system evolve before measuring the observ- 
ables) there is no evidence of a divergence of the spin- 
glass susceptibility with the size of the systenS. This 
is at variance with the results obtained in case of short- 
range Edwards- Anderson spin glasses, where the grow-iag 
of the correlation length has been generally observed 2 ^. 

IV. OFF-EQUILIBRIUM DYNAMICS 

In the previous section we have seen that the relaxation 
time increases very fast as the temperature is lowered, 
and it is equally well fitted by a finite-time Vogel-Fulcher 
law and by the law t — A exp(B /T 2 ) . Even in this second 
hypothesis, in the whole low temperature range T < 1, 
the time scales we can computationally reach for reason- 
ably large systems are by far shorter than the equilibra- 
tion time. History dependent effects and aging are then 
to be expectedcil 

Recent results in the study of the dynamics following a 
quench in the low temperature phase in spin glasses, give 
prediction about the scaling of the correlation function 
and the linear response function for large times in the 
off-equilibrium regime, and on the dependence of these 
quantities on the. time t w that the system has spent at 
low temperature!! These functions are defined respec- 
tively as 

C{t, t w ) = i 5^(o-i(f + t w )ai{t w )) 

where hi is a local magnetic field applied to the system. 

For a complete exposition of the off-equilibrium theory 
of the glassy dynamics we refer the reader tcotl Here we 
limit to reassume briefly some features relevant to the 
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present discussion. For large t w the following scenario is 
found: there is a first regime, for small t (t << r(t w ) 
see below) where the dynamics has features similar to 
those of an equilibrium system. The correlation func- 
tion is independent of t w in this regime, and the re- 
sponse function is related to the correlation function by 
the fluctuation-dissipation theorem. In this regime the 
correlation function monotonically decreases from 1 to 
a value qEA, which defines an off-equilibrium parameter 
analogous of the Edwards- Anderson parameter. In addi- 
tion to this equilibrium-like regime, in the class of models 
of interest for this paper, there is a regime in which the 
correlation decays from qEA to zero, and the correlation 
has the scaling form 

C(t,t w ) = C{t/r(t w )). (9) 

The "effective relaxation time" r(t w ) is an increasing 
(and diverging) function of t w , that the theory is not 
able to predict, and seems to be rather si^tem depen- 
dent. In some cases it is found r(t w ) = t w c3c3 : we will 
see that this is not the case for the present model. 

As far as the behaviour of the response function is con- 
cerned, it is found that the function 

r. . \ _ TR(t, t w ) 
x\h i-w) — gc(t,t w ) ac(t,t w ) ' ^ V) 

dt w dt 

sometimes called fluctuation-dissipation ratio, depends 
on its time arguments in a quite special way: through 
the dependence on t and t w of the correlation function 
itself, 

x(t,t w )=x(C(t,t w )). (11) 

The value x = 1 corresponds to the Fluctuation- 
Dissipation Theorem relation and it is valid for qEA < 
C(t,t w ) < 1. A non zero value of x(C) is found when 
aging effects are present in the response function, and 
testifies the memory of the system about its history. In 
the mean- field p-spin model it is found that x{C) is equal 
to a constant x smaller than 1 in the whole interval 
< C < qEA- For a discussion of the behaviour of 
this quantity in the Edwards- Anderson model, as well as 
for a qualitative discussion of its behaviour in a different 
glassy scenario seeEJ. 

The response function is measured from simulations of 
"zero-field cooled experiments"E3. Starting at time zero 
from a random configuration, we let evolve the system at 
constant temperature in zero field for a time t w . At t w 
we switch on a small magnetic field h, and measure the 
relaxation of the magnetization as a function of time. In 
the linear response regime, the magnetization m(t,t w ) is 
given by 

rt+t w 

m(t,t w ) =h J ds R(t,s), (12) 

a relation that, assuming the validity of (^lj), takes the 
form 



h /* ^ h 
m{t, t w ) = - dq x{q) = - X (C(t, t w )). (13) 

1 JC(t,t w ) 1 

The mentioned behaviour of x(C) reflects in 

X{C) = { 1 - q EA (l -x)-Cx °C <q EA (W) 

It has been noted i ni that, while a scaling behavior of 
the kind (^|) is common to a glassy behavior, and phenom- 
ena of domain-growth in phase separations, the function 
x(C) seems to be non zero only in the aging regime of the 
glassy systems. In order to discriminate glassy behavior 
from domain-growth like mechanisms another quantity 
has been recently proposecO. This is defined consider- 
ing the evolution of two replicas of the system, {(Ji\ and 
{t;} which follow identical evolution up to a time t w , and 
independent evolutions afterwards. One then considers 
the correlation 

Q{t,t w ) = (ai(t + t w )Ti(t + t w )), (15) 

that, by definition, is equal to one for — t w < t < 
0. Barrat, Burioni and Mezard have discussed in de- 
tail the meaning of this variable proving that g — 
limi^^oo lim^oo Q(t, t w ) is different from zero in some 
domain-growth models, while it tends to zero in mean- 
field spin-glasses. This shows that, in the last case, 
the two typical trajectories explore different regions of 
the phase space. In equilibrium, it holds the relation 
Q(t/2,t w ) — C(t,t w ) (with both quantities independent 
of t w ). In has been also shown that the same relation 
holds in a trap model even in the aging situation. 

We have investigated the behavior of the functions C, 
R and x i n numerical simulations for a large size L = 30 
and temperatures T = 0.5,0.7. The two temperatures 
give qualitatively similar results. In figure 7 we plot our 
data for C(t,t w ) and Q(t,t w ). As anticipated, the scal- 
ing r(t w ) = t w is not obeyed. A best fit of the form 
r(t w ) — gives a = 0.77 and produce a fairly good col- 
lapse of the data on the same master curve. We did not 
investigate a possible dependence of a on the tempera- 
ture. The behavior of Q(t, t w ) indicates that the parame- 
ter g is zero in this model as predicted by the mean-field 
theory. We see that the relation Q(t/2,t w ) = C(t,t w ) 
is obeyed at short times, and does not work so badly 
even at large times. However a rescaling of the kind 
Q(</(1.5), t w ) — C(t,t w ) fits better the data. 

As far as the function x(C) is concerned, we plot in 
figure 8 the rescaled magnetization \ = (T/h)m(t,t w ) 
versus C(t,t w ) for different t w . The figure shows clearly 
that for the waiting times here considered we are vary far 
from an asymptotic regime where \ is independent of t w . 
However, the slope of the curve for small C seems not to 
vary too much with t w , and it is roughly equal to x = 0.4. 
The region where x(C) = 1 — C terminates roughly at a 
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value qEA = 0.97. Assuming (as supported by mean-field 
theory) that the Edwards-Anderson parameter defined 
in this way and the one defined in statics take equal or 
comparable values, we argue that the quasi-states domi- 
nating the thermodynamics are quite "narrow" , and, as 
announced, the intra-state entropy is quite small. To 
our knowledge, the function x(C), or equivalently x(C) 
has never been measured in experiments. Its determina- 
tion, which would involve independent measures of a time 
correlation function and its associated response function, 
would be a good test of the spin-glass scenario in glasses. 



V. CONCLUSIONS 

In this work we have investigated the glassy behav- 
ior of a short-range disordered p-spin interaction model. 
We have focused to the case p — 4 in three dimen- 
sions, reporting for comparison purpose some results for 
the two-dimensional system. The mean-field version of 
this model displays a static transition where the config- 
urational entropy is nearly zero and replica symmetry 
breaks. We have found that the finite dimensional model 
shows interesting features, archetypical of real labora- 
tory glasses. Both in two and in three dimensions the 
equilibrium autocorrelation function follows a stretched 
exponential form at low enough temperature. While in 
two dimensions the relaxation time increases according to 
the Arrhenius law, in dimension three our data demon- 
strate a faster increase. This is compatible both with 
a Vogel-Fulcher law, and with an exp(A/T 2 ) behavior. 
In favor of the finite-temperature transition scenario, we 
can only offer the extrapolations of the high temperature 
data of the entropy. A better support to the theory is 
furnished by the relation t ~ exp(A/(TS(T)) that we 
observe in the temperature window we explored. The 
search in several quantity, of a growing statical corre- 
lation length associated to the growing relaxation time 
gave negative answer. In all cases we found very small, 
and nearly temperature independent correlation lengths. 
Of course we can not exclude that the correlation length 
begins to grow at temperatures lower that the ones at 
which we could equilibrate the system. 

The study of the off-equilibrium dynamics of the 
system confirmed the qualitative features predicted by 
Mean-Field Theory. The behavior of the auto-correlation 
function shows aging. As times goes by the dynamics 
becomes slower. The effective relaxation time, defined 
as the characteristic time for the correlation function to 
vanish, grows as a power of the waiting time. The anal- 
ysis of the function Q(t,t w ) shows that typical trajecto- 
ries which coincide at time t w explore uncorrelated re- 
gions of the space at later times. The behavior of the 
fluctuation-dissipation ratio, displays the qualitative fea- 
tures expected, showing long-range memory effects in the 
aging regime. However the residual dependence on t w in- 
dicates that at the times we have reached the eventual 



asymptotic behavior is still very farcj. 

The evolution of glassy systems is often described phe- 
nomenologically as rare jumps in a landscape of "traps" 
(metasiable states) inside which most of the time is 
spent£3. This picture matches quite well with the theo- 
retical idea of the metastable states destabilized by acti- 
vated processes. It has been pointed out inE3 that jumps 
among uncorrelated traps imply the relation Q(t,t w ) = 
C(2t,t w ). Our data show Q(t,t w ) is indeed close to 
C(2t,t w ). If this would be confirmed by more precise 
and systematic studies, would be an evidence in favor of 
the trap models. In this respect it would also be useful 
the analysis of the self-averaging properties of some local 
quantity. 

We would like to conclude mentioning some of the 
problems left open by this work. For example it would 
be interesting to study the behavior of the dynamic 
Edwards- Anderson parameter as a function of tempera- 
ture, or the large-time decay of the energy to its asymp- 
totic value. We have also seen, that the equilibrium re- 
laxation in high temperature does not seem to proceed in 
two steps (/3 and a) as it commonly observed in labora- 
tory, and as it is predicted by Mode Coupling Theory. As 
far as we know, the /3-relaxation process has never been 
observed in a spin system. We do not know if this is just 
due to a difficulties to disentangle the two processes due 
the high value of the Edwards- Anderson parameter in the 
neighborhood of the transition, or to the real absence of 
the (3 process. We think however that the understanding 
of this point can shed some light on the nature of the 
/^-relaxation process. On the theoretical side, it would 
be very interesting to understand if the observed value of 
the exponent 7 = 1 in three dimensions, in contrast with 
the replica theory prediction for the Potts-glass 7 = 2 
has to be imputed to the plane inversion symmetry, or if 
it is observed even in absence of this. In that respect sim- 
ulations of the Potts-glassEl could give important hints. 
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Figure Captions 

Fig. 1. The energy as a function of the temperature for 
different cooling rates r. It is apparent the calori- 
metric glass transition around T w 0.7. The region 
around the transition is magnified in the inset. 

Fig. 2. The high-temperature entropy (data was collected 
for one sample with lattice size L = 10) versus the 
inverse temperature, together with the rational fits 
mentioned in the test. The fitting parameters are, 
respectively for the s\ and S2 forms: a = —0.302, 
b = 1.866 and a = -0.039, b 2 = 2.060, b 4 = 1.550. 

Fig. 3. The correlation function for various temperatures: 
T = 2.8 (□), T = 2.4 (x), T = 2.0 (A), T = 1.6 
(0)> T = 1.2 (+). The lines are stretched expo- 
nential fits. 

Fig. 4. The relaxation time as a function of the temper- 
ature toghcther with the Vogel-Fulcher fit (dotted 



line) and the exp(/3 2 ) fit (slashes). The inset shows 
the temperature dependence of the exponent b. 

Fig. 5. The inverse of the logarithm of the relaxation time 
versus the TS{T). The linear dependence is the 
one predicted by the Adam-Gibbs form. 

Fig. 6. Spin-glass susceptibility \ f° r different sizes L = 
4, 5, 6, 7 averaged over 100 samples versus the tem- 
perature. There is no evidence of a growing corre- 
lation length. 

Fig. 7. The functions C and Q at T = 0.5. We plot on 
the same graph C(t,t w ) and Q(t/2,t w ). We see 
that the relation C(t,t w ) = Q(t/2,t w ) is quite well 
obeid even for small values of C. The curves for 
t w = 10, 100, 1000, 10000 are depicted. 

Fig. 8. The function x(C) for T = 0.7 and waiting times 
t w = 10, 100, 1000. The curves show a strong de- 
pendence on t w indicating that we are far from the 
eventual asymptotic regime. 
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